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We study the electron-phonon interaction in the strongly correlated superconducting cuprates. 
Two types of the electron-phonon interactions are introduced in the t — J model; the diagonal 
and off-diagonal interactions which modify the formation energy of the Zhang-Rice singlet and its 
transfer integral, respectively. The characteristic phonon-momentum (q) and electron-momentum 
(k) dependence resulted from the off-diagonal coupling can explain a variety of experiments. The 
vertex correction for the electron-phonon interaction is formulated in the SU(2) slave-boson theory 
by taking into account the collective modes in the superconducting ground states. It is shown that 
the vertex correction enhances the attractive potential for the d-wave paring mediated by phonon 
with q — — S), 0) around 5 = 0.3 which corresponds to the half-breathing mode of the oxygen 
motion. 

PACS numbers: 74.20.-z, 74.25.Kc, 74.20.Rp, 74.25.Ha 



I. INTRODCUTION 

It is widely accepted that the central issue in the high 
transition-temperature (Tc) superconducting (HTSC) 
cuprates is physics of the doped Mott insulator. The 
Mott-Hubbard (charge-transfer) insulating state with the 
antiferromagnetic long-range order in the parent com- 
pounds is well understood in strong Coulomb interaction 
between electrons. It is, then, well recognized that the 
electronic models with strong electron correlation provide 
a good starting point to reveal the electronic structure in 
doped Mott insulator. For example, one of the successful 
examples is the slave-boson mean-field theory in the two- 
dimensional t — J model and its extensions to the gauge 
theory which can explain a wide range of unusual proper- 
ties both in the normal and superconducting states; the 
phase diagram in hole concentration x and temperature 
T, paring symmetry, pseudogap and so on. 15 

In spite of such successful results in the correlated elec- 
tronic models, one shall have the questions that the elec- 
tronic model alone is enough to explain essential physics 
in HTSC cuprates. Actually, since the early stage of 
the HTSC researches, a variety of experiments exhibit 
the lattice/phonon anomalies and the strong electron- 
lattice coupling which stimulate a number of theoreti- 
cal works. 6 " 11 The finite isotope effects on Tc is one of 
the direct evidence of the phonon contribution to the 
paring interaction. The oxygen isotope coefficient ao 
in La2- x Sr :E Cu04 (LSCO) around x = 1/8 exceeds 0.5 
being greater than the expected value in the Bardeen- 
Cooper-Schrieffer (BCS) theory. 12 The isotope shift is 
also found in the magnetic penetration depth at zero 
temperature A(T = 0) (Ref. 13) which can not be un- 
derstood within the standard BCS scenario. Anomalous 
phonon dispersion relation is directly observed by the 



inelastic neutron scattering experiments. A strong soft- 
ening and broadening of the highest longitudinal optic 
(LO) phonon with the oxygen in-plane Cu-0 bond vi- 
brations occurs along the (q x ,0) direction in LSCO. 14:15 
By doping of holes, the phonon energy (~ 85meV) dras- 
tically reduces about 20% in the zone boundary, and is 
recovered in overdoped compound. 16 

A recently observed kink structure in high-resolution 
angular resolved photoemission spectroscopy (ARPES) 
spectra 17-20 triggers reexamination of the strong 
quasiparticle-phonon interaction in cuprates. 21-25 Here 
we summarize characteristic nature of this kink structure: 
(a) The kink structure is experimentally confirmed in a 
variety of the p-type cuprates: LSCO, Bi 2 Sr 2 CaCu208 
(Bi2212), Bi 2 Sr 2 CuO 6 (Bi2201), (Ca 2 _ ;c Na a: )Cu0 2 Cl 2 
(Na-CCOC), Tl 2 Ba 2 Cu0 6 (T12201) and so on. On 
the other hand, this structure is weak in the n-type 
Nd 2 _ x Ce x Cu0 4 (NCCO). (b) The kink energy uj kmk 
about 70meV is almost universal. The quasi-particle ve- 
locity ratio above and below LOkink, he. R v = v(lu > 
<^kink)/v(u> < Wkink), is about 2 in optimal doped p- 
type cuprates and does not show the strong momentum 
dependence along the fermi surface. Below u> k ink, the 
quasi-particle width decreases rapidly. This is consis- 
tent with the rapid drop of the quasi-particle scattering 
rate deduced from the optical conductivity data. 26,27 (c) 
The kink structure is seen far above Tc, and remark- 
able changes across Tc are not observed, (d) The quasi- 
particle velocity below the kink v(u> < u>ki n k) along the 
nodal direction is almost universal in a wide range of 
HTSC. 28 

This phenomena indicate the strong coupling between 
a bosonic excitation with energy Wkink and the quasi- 
particles. It is suggested that there is a clear correla- 
tion between the quasiparticle velocity ratio across the 
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kink R Vl related to the electron-boson coupling strength, 
and the superconducting gap amplitude. 21 A scenario 
based on the magnetic-resonance mode around 41meV 
as a bosonic excitation 19 ' 29 ' 30 may be at a disadvan- 
tage; this excitation is only confirmed below Tc in 
YBa 2 Cu 3 6+;c (YBCO), Bi2212 and T12201. One of the 
plausible candidates for the origin of the kink structure 
is the optic phonon. The kink energy uJkink is close to 
that of a LO phonon with the oxygen in-plane vibration 
whose anomalous softening and broadening are observed 
as described previously. This phonon scenario is strongly 
supported by the recent ARPES experiments for the iso- 
tope effects on the kink; in a Bi2212 sample with the 
O 18 isotope, the kink energy u>kink decreases about few 
meV from the standard O 16 results. This result corre- 
sponds to the simultaneously observed lower shift of the 
LO phonon Raman spectra in the O 18 sample. 31 

Here wc mention another kink structure in the ARPES 
spectra in Bi2212. 32-34 This is found near (tt, 0) and 
(0, 7r) and, its kink energy is about 40meV measured from 
the Fermi energy. It has been claimed that this kink 
structure almost disappears above Tc, but recent more 
extensive study 34 found that it persists even above Tc. 
Therefore the c-axis buckling phonon, which has the ap- 
propriate energy of 40meV, is the most promising among 
the various candidates such as the magnetic resonance 
mode and the superconducting gap. However, we focus 
below the half-breathing mode which plays the major role 
in the nodal region. 

In this paper, we study the electron-phonon interaction 
in the strongly correlated HTSC cuprates. The interac- 
tions between electron and phonon of the in-plane oxygen 
vibration mode are formulated in the t — J model. There 
are two types of the couplings, i.e. the off-diagonal and 
diagonal interactions which modify the inter-site hopping 
of the Zhang-Rice singlet and its formation energy, re- 
spectively. We calculate the effective paring interaction, 
the quasiparticle rcnormalization factor, the ARPES, 
tunneling and optical spectra for the off-diagonal and 
diagonal couplings where the fermionic degree of free- 
dom is assumed as a quasiparticle with the renormalized 
band energy. We show that the characteristic phonon- 
momentum (q) and electron-momentum (k) dependences 
resulted from the off-diagonal coupling are quite consis- 
tent with the experiments. This formulation and numer- 
ical results are presented in Sec. II. The vertex correction 
in the electron-phonon coupling is known to play a cru- 
cial role in correlated systems. 40 ' 41 This is studied, in the 
slave boson picture, by taking into account the fluctua- 
tions around the mean-field saddle-point solutions, that 
is, the d-wave paring order parameter A, the fermion 
hopping x, the Lagrangian multipliers a and the bosonic 
field h. We formulate this vertex correction in the SU(2) 
slave-boson theory 35-38 which respects an exact SU(2) 
gauge symmetry in the undoped case, and is expected 
to provide a better starting point than the U(l) theory 
in the underdoped region. The fluctuations are treated 



as the collective modes in the ordered states, 39 and the 
electron-phonon interactions are corrected by the interac- 
tion between the collective modes and fermions. We show 
that the vertex correction enhances the effective interac- 
tion for the d-wave paring around q = (ir(l — 8), 0) with 
8 = 0.3 where the remarkable softening and broadening 
of the LO phonon are observed in the inelastic neutron 
scattering experiments. The roles of the vertex correc- 
tions are examined in Sec. III. Section IV is devoted to 
the summary and discussion. A part of the present the- 
oretical works were briefly discussed in Rcf. 21. 



II. ELECTRON PHONON INTERACTION IN T-J 
MODEL 

A. Formulation of the Electron-Phonon Interaction 

We will consider the electron-phonon interaction in the 
t - J model: 42 
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where Ci a is the annihilation operator of a hole at site i 
with spin cr, and is defined in the Hilbert space excluding 
the double occupancy. Si is the 5=1/2 spin operator at 
site i. The formation energy of the Zhang- Rice singlet 
at site i and its transfer integral Uj between i and j are 
represented by the energy parameters in the original p — d 
model as 
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and 
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respectively. The charge-transfer energy A pc i(i) is defined 
to be the energy difference of the p and d levels at site 
i: A p d(i) = s p — Sd(i) in the hole picture, and t p d is the 
transfer integral between the nearest neighboring (NN) p 
and d orbitals. We assume that U — A p d(i) >> A p d(i). 

Here, we introduce a motion of the oxygen ions along 
the Cu-0 bond in the Cu02 plane corresponding to the 
high-frequency optical phonon of our interest. The dis- 
placement of the negatively charged O ion along the Cu- 
O bond changes the electro-static Madelung potential 
acting on a Cu site and shifts the energy of the d level as 
£d(i) = £do + Se d (i) with 



5£d{i) = gu(i), 



(4) 



within the linear electron-lattice coupling. u(i) is a linear 
combination of the O displacements around the Cu site 
i (Fig. 1(a)) defined by 
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u(i) =u x (i + y) - u x (i - y) 
+ u y (i + y) - u y (i - y) ; 



(5) 



where Ui(i ± ^f) (I = x,y) indicates the displacement 
of the ion at site i + tt along the direction I. The 
coupling constant g is positive, because the approach of 
the negatively charged O ion to the Cu site lowers the 
energy level for a hole. The O p level is not changed 
within the linear electron-lattice coupling, i.e. Se p = 0, 
because of symmetry. The modification of the charge- 
transfer energy is then given by A pc i(i) = Ao — gu(i). 
The O ion displacement also changes the transfer integral 
between p and d orbitals as 



t pd = t ± g t ui (i ± ^) 



(6) 



for a bond connecting a Cu site i and an O site i ± ai/2. 
The coupling constant gt has an opposite sign of to. As 
a result, the modulation of the formation energy of the 
Zhang-Rice singlet is given by 
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which is termed the diagonal electron-lattice coupling. 
We note that the two terms in the right-hand side have 
opposite sign with each other. In the similar way, the 
transfer integral of the Zhang-Rice singlet is modified as 
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tij = -f L +9off{u(i) + u(j)}, 
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termed the off-diagonal electron-lattice coupling. It is 
worth noting that the modulation of t p d does not change 
tij. This is because, with a shift of a O ion at site i + ai/2 
along the positive I direction, a decrease in t p d between 
site i and site i + ai/2 is canceled out in the linear order 
by an increase in t p d between i + a; and i + ai/2. The 
electron-phonon interaction in the t — J model is also 
brought about by the modulation of the superexchange 
interaction Jjj due to the O displacements. This type of 
electron-phonon coupling was first suggested to explain 
the two-magnon absorption spectra in undoped cuprates 
in Ref. 43. Here, we consider this electron-phonon cou- 
pling in the framework of the slave-boson mean-field the- 
ory. The detailed formulation based on the salve-boson 
theory is presented in Sec. Ill A. The constrained opera- 
tor Ci a is represented, in this picture, by a product of the 
spinon (fermion) / i(T and holon (boson) hi as 



(a) 



u [i+a /l) 



O O / O O j 

— u x (i-a x /2) u x (i + a t /2) 

-u {i-a /2) 



9 



(b) 
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FIG. 1. (a) A linear combination of the O ion distortion 
u(i) = u x (i + a x /2)+u y (i + a y /2) —u x (i — a x /2) — u y (i — a y /2). 
Open and filled curcles indicate the O and Cu ions, respec- 
tively, (b) Schematic view of the half-breathing mode of the 
O vibration. +, — and indicate signs of the modulation 
of the transfer integral Uj between the NN Cu ions in this 
displacement. 
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with the local constraint ^2 a fj a fia + h\hi = 1 in or- 
der to exclude the doubly occupied states. In the mean- 
field approximation, the J term in the t — J model is 
rewritten by using the fermion hopping order-parameter 

Xij(= Ysaftfjer) as 



J ij 



Si 



JijXijXij 
J, 



^-^((xljhij+xUxii)), (12) 



which is a similar form to the t term in Eq. (1). The 
fcrmion-lattice coupling through the modification of J is 
then reduced to the off-diagonal coupling as 



Jij = Jo + gj {u(i) + u(j)} , 
where the coupling constant is given by 



9J 



(13) 



(14) 



In the standard slave-boson mean-field scheme, where 
the boson amplitude is taken to be \fx in the ground 
state, the diagonal coupling constant gdia and the off- 
diagonal one from the modulation of the transfer integral 
g ff are scaled by the factor x, and the coupling caused 
by the superexchange interaction gj has a factor 1 — x. 



3 



This implies that, in the Mott insulating limit of x — > 0, 
due to suppression of the charge fluctuation, the diag- 
onal coupling becomes irrelevant, while the off-diagonal 
type survives because the modulation of J affects the 
spin channel. Even at nonzero (but small) x where the 
diagonal coupling becomes relevant, the vertex correc- 
tion for the electron-phonon interaction, induced by the 
holon fluctuation, suppresses the diagonal coupling at the 
large momentum transfer. 40 ' 41 The vertex correction for 
the off-diagonal coupling is presented in Sec. Ill in more 
detail. 

Consequently, the electron-phonon interaction Hamil- 
tonian in the t — J model is summarized as 



k,q,cr P 



fl f, (bZ + Vl) 
J k+qcr J k<y \ q T -q J 



(15) 



where b*i is the annihilation operator of phonon with mo- 
mentum q and the mode \i. We take the two independent 
modes [i — x and y where the O ions vibrate along the x 
and y directions, respectively. The electron-phonon cou- 
pling constant is given by a sum of the two-kind couplings 



9^k,q)=9" dia (k,q)+9 , : ff (k,q), 



(16) 



where 
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+ cos(k x + q x ) + cos (k y + %)}, 



(18) 



with the oxygen mass M and the phonon frequency ujt. 
The off-diagonal coupling constant g f / is redefined that 
g ff includes the contribution from g.j in Eq. (14). The 
diagonal-coupling constant gdia(k,q) docs not depend 
on the electron momentum k as usual. Here, we note 
the characteristic electron-momentum k and phonon- 
momentum q dependences of the off-diagonal coupling 
goff(k,q): (i) the fcrmion degree of freedom does not 
couple to the two-dimensional oxygen breathing mode, 
i.e. ggff{k,q = (7r, 7r)) = 0. This is attributed to the 
fact that, for the breathing mode, u(i) + u(j) in Eq. (9) 
becomes zero in all the NN Cu-Cu bonds, (ii) The cou- 
pling constant for the half-breathing mode where the O 
ions shift along / does not depend on the I component of 
the electron-momentum, e.g., g^Ak,q — (n, 0)) is inde- 
pendent of k x . This is because this O ion displacement 




(0, 0) k (it, 0) 

FIG. 2. Contour map of the effective paring interaction 
F(q) in the off-diagonal coupling case (a), and in the di- 
agonal coupling case (b). The energy gap is chosen to be 
4>o9off = 0.035eV in (a), and (polldia = 0.035eV in (b). Nu- 
merical data are plotted in a unit of eV. The phonon energy 
is chosen as w p / l =0.07eV. 



does not change Uj along the direction I but change Uj 
perpendicular to I. This characteristics of g„ff(k, q) arc 
schematically shown in Fig. 1 (b). 



B. Physical consequences of the off-diagonal and 
diagonal electron-phonon couplings 

We now turn to the numerical results for various phys- 
ical quantities calculated from the electron-phonon inter- 
action Hamiltonian introduced above. We pay attention, 
in particular, to the implications of the characteristic mo- 
mentum dependence of the off-diagonal electron-phonon 
coupling. Then, we treat, in this section, an electron as a 
quasiparticle described by the renormalized energy band 
£r. This treatment corresponds to the mean- field approx- 
imation in the slave-boson theory where the Lagrangian 
is given in Eq. (35) in Sec. Ill A. In spite of the detailed 
form for the quasiparticle energy £g in the slave-boson 
theory (see Eq. (36)), £g is determined by the tight-bind 
fitting in Bi2212. The six tight-binding parameters are 
adopted by following the results in Refs. 29 and 44. The 
fluctuations from the mean-field saddle-points leading to 
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the vertex correction for the electron-phonon interaction 
are studied in the next section. The energy parameter 
values used in this section are given in a unit of eV. We 
assume that the phonon frequency Lu'i is independent of 
the mode /i and the momentum q as to- — u> p h, and 
choose uoph =0.07cV. The electron-phonon coupling con- 
stants are scaled as g d ia(off) = 9dia(off)l ^IMNut, and 
the damping constant is introduced in the calculation ol 
spectra as 7 — 0.005eV. 

We first show the effective paring interaction for the 
d-wave superconductivity due to the electron-phonon in- 
teraction. The momentum dependence of the effective 
electron-electron interaction is derived by integrating out 
the phonon degree of freedom in the electron-phonon in- 
teraction Hamiltonian. By using the second-order per- 
turbational processes with respect to the electron-phonon 
coupling, the effective Hamiltonian is obtained as 

H *ff = -^ E E 9^Q) 9 ^k',-q) 
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In the mean-field approximation according to the BCS 
theory, we obtain 



FIG. 3. Contour map of the fermion spectral density 
A(k,u>) for the off-diagonal case (a) and (b) (<? o //=0.03eV), 
and for the diagonal case (c) and (d) (gdia~0.02oV). Solid 
lines indicate the spectra with Ao = and dotted lines in 
(a) and (b) are for the spectra with Ao = 0.035eV. Gray 
bold lines indicate the bare quasiparticle dispersion £g. The 
ft f- ft fr, ,. (19) phonon energy is uj ph =0.07eV. Numerical data are plotted in 

' k+qa' Ka k' —qa' K 17 v ' „ , ,_1 
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unit of eV 

strongly suppresses the d-wave paring. That is, this di- 
agonal coupling competes with the magnetic paring in- 
teraction. 

The ARPES, tunneling and optical spectra are ob- 
tained in the Nambu-Eliashberg formulae. We calculate 
the self-energy for the 2x2 fermion Green's function: 



9,M 



H(k,iuj n ) 



where |£g — ^fc+gl <<; ^ 1S assumed. The paring or- 
der parameter is introduced as = (f^f^) and is 
parameterized as 0£ = (j) (cosk x — cosfc y )/2 for the d- 
wave paring. F(q) indicates the effective paring inter- 
action between electrons contributed from the phonon 
with momentum q. In Fig. 2(a), the contour map of 
F(q) for the off-diagonal coupling is presented. The neg- 
ative (positive) region corresponds to the attractive (re- 
pulsive) interaction for the d-wave paring. It is found 
that the phonon of the momentum (q x ,0) around q x = n 
strongly contributes to the d-wave paring. This phonon 
corresponds to the half-breathing mode of the O vibra- 
tion where the anomalous softening of the dispersion re- 
lation is observed in the neutron experiments. 14,15 The 
vanishing F(q = (7r,7r)) indicates that the off-diagonal 
electron-lattice coupling does not disturb the attractive 
magnetic interaction dominated around (ir, n). This orig- 
inates from the characteristic off-diagonal coupling that 
the breathing-mode of the O vibration does not couple to 
the fermion as pointed out above. The paring function 
F(q) for the diagonal coupling is presented in Fig. 2(b) 
where the attractive interaction appears in the region of 
Qx + Qy ^ I"- The most remarkable difference from the 
off-diagonal results is seen that the phonon around (ir, n) 



^2 T ^G {k - q, iu n - iiv m )T 3 



x\g»(k,q)\ 2 DZ(q,iu m ), (21) 
where we introduce the bare fermion Green's function 

iLU n T + £gT3 + A^Tl 



G (k,iui n ) 
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with 
tion 



Dg(q,iw m ) 



£i + A- and the bare phonon Green's func- 

' k k 



lUJ r 



IT] 



(23) 



t; (I = 1,2,3) is the Pauli matrices, t° is a unit ma- 
trix, and Ar(= Ao(cosfcz — cosk y )/2) is the supercon- 
ducting energy gap. The energy and momentum depen- 
dences of the occupied spectral density corresponding to 
the ARPES spectra are obtained by the (l,l)-component 
of the Green's function: A(k,ui) = — ^lmGu{k, iu> n — > 
uo + iS). The calculated results are presented in Figs. 3 
(a) and (b) for the off-diagonal coupling case, and in 
(c) and (d) for the diagonal case. The ARPES spec- 
tra show a kink structure below which both the spectral 
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FIG. 4. Fermion density of states N(lo). The diagonal and 
off-diagonal coupling constants (ffo//, <?t«a) are chosen to be 
(0.03eV, 0) for a solid line, (0.015eV, O.OleV) for a broken 
line, and (0, 0.02eV) for a dotted line. Other parameters are 
uj ph = 0.07eV and A = 0.035eV. 



intensity and the quasi-particle lifetime remarkably in- 
crease. The kink energy in the superconducting state is 
less than A + cj p h because of the d-wave gap. Although, 
at a glance, the kink effect in the spectra seems to be 
remarkable in the anti-nodal direction, this impression 
comes from the almost flat dispersion around (n, 0) of 
the bare energy band. As shown later, the quasiparticle 
rcnormalization factor Z(k,Lu) dominating the velocity 
change at the kink energy is slightly weaker in the anti- 
nodal direction for the off-diagonal coupling case. This 
is consistent with the recent ARPES experiments where 
the quasiparticle velocity does not show a huge change 
at LOkink around the anti-nodal direction, unlike the pre- 
diction in the magnetic mode calculations. 20 ' 29 
The tunneling density of states 



N(uj) 



^ Tr ImG(fc, iu>„ — ► u + iS), 



(24) 



are presented in Fig. 4. An asymmetrical shape of the 
spectra is attributed to the almost flat band around the 
(it, 0) point below the fermi energy. Outside the super- 
conducting gap, a dip and hump structure appears, be- 
ing consistent with experiments. 45 ' 46 In the off-diagonal 
coupling case, the weak hump structure is located at 
±(A +Lo p h)- With increasing the diagonal-coupling pa- 
rameter, this structure becomes pronounced, in particu- 
lar, in the negative side, and the position of the structure 
shifts to the lower energy. It is seen simultaneously that 
for the strong diagonal coupling, the gap structure be- 
comes shallow and is gradually collapsed. 

The quasiparticle scattering rate deduced from the op- 
tical conductivity is calculated in the memory-function 
formalism. The optical conductivity is represented by 
utilizing the memory function M aa (uj) which corresponds 
to the inverse of the energy-depend life-time r(u) as 
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FIG. 5. Imaginary part of the memory function ImM(cj). 
The diagonal and off-diagonal coupling constants (g ff,^dia) 
are chosen to be (0.03eV, 0) for a solid line, (0.015eV, O.OleV) 
for a broken line, and (0, 0.02eV) for a dotted line. Other 
parameters are lu v h = 0.07eV, Ao=0 and t = 0.3eV. We use 
that xo = 0.81ta 2 calculated for the half-filled tight-binding 
model in the two-dimensional square lattice with NN electron 
hopping. 



a aa (uj) 



47T w + M aa (u) ' 
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with the plasma frequency u>p and the Cartesian coordi- 
nate a. The memory function M aa (u) is defined by 47 

M aa (L0) = — { Xqq M -Xaa(0)}. (26) 

Xao(w) is the Fourier transform of the retarded Green's 
function of the operator A a defined by the equation of 
motion for the current operator j a as A a = [j a , H] . Ex- 
plicitly, 

A a — — 2ta ^ {sina(/c Q + q a ) — sinafc a } 



(27) 



where we consider the current between the NN sites. In 
the calculation of M aa (oj), the correlation functions of 
fermion and phonon operators are evaluated for the non- 
interacting Hamiltonian. xo i n Eq. (26) is the static limit 
of the response function of j a . The calculated ImM(w) 
are shown in Fig. 5. lmM(uj)(=lmM xx (uj) =lmM yy (u))) 
is rapidly depressed below Lo p h- Introducing the diag- 
onal coupling, the scattering rate increases above the 
threshold energy and its reduction below the energy be- 
comes steep. This large depression in the quasiparticle 
scattering rate deduced from the optical spectra is com- 
monly observed in a variety of HTSC, e.g. LSCO, YBCO, 
T12201 and NCCO around w = 400 - 800cm- 1 (Ref. 27). 
This threshold energy is rather larger than the so-called 
pseudogap estimated in the inelastic neutron scattering, 
ARPES and tunneling spectra, and is close to the kink 
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energy related to w p h in the present scenario. Absolute 
values of the calculated ImM(u) are compared with the 
experimental quasiparticle scattering rate, xo is esti- 
mated in the two-dimensional tight-binding model with 
the NN electron hopping. Around the threshold energy, 
the experimental t(oj) -1 is about 2000cm -1 in under- 
doped Bi2212 and LSCO with x = 0.14, 1000cm- 1 in 
NCCO, and 500cm- 1 in overdoped Pb doped Bi2212. 27 
These values are comparable to the present results of 
ImM(cj) in Fig. 5. That is, the electron-phonon interac- 
tion provides a large portion of the scattering rate in the 
optical region. 

The memory function ImM(w) at u> = corresponds to 
the scattering rate in the DC reisistivity. The tempera- 
ture dependence of ImM(o> = 0) is presented in Fig. 6. It 
is known that the scattering rate deduced from the exper- 
imental resistivity data is of the order of 2T. 48 The cal- 
culated ImM(0) gradually increases around w p h/5 with 
increasing temperature, and its absolute value is rather 
smaller than 2T. These results indicate that the phonon 
contribution is less dominant for r(O)- 1 in the resistivity, 
in contrast to the scattering rate in the optical region as 
mentioned above. It is noted that the superconducting 
gap Ao is put to be zero in Figs. 5 and 6. When we 
introduce the finite Ao to describe the pseudo-gap state, 
ImM(0) will be further reduced due to the reduced final 
density of state. 

The momentum dependences of the fermion renormal- 
ization factor Z(k,co) = (1 — <9E(fc, ivj/div)^ 1 at ui = A 
are presented in Figs. 7(a) and (b) for the off-diagonal 
and diagonal couplings, respectively. The broken lines 
indicate the fermi surface for £g. In the off-diagonal cou- 
pling case, ReZ(k,ui = Ao) is close to one around the 
k x + k y = 7r line in the Brillouin zone, and is strongly re- 
duced around (0,0) and (tt,tt). In contrast, the momen- 
tum dependence is much weaker in the diagonal case, 
reflecting that the coupling constant g^ ia (k,q) is inde- 
pendent of the electron-momentum k as seen in Eq. (17). 
The calculated results for the off-diagonal coupling can 
explain the anisotropic quasiparticle scattering rate sug- 
gested in the recent ARPES experiments: Full width 
half maximum (FWHM) of the quasiparticle energy dis- 
tribution curve is measured along the fermi surface in 
Pb-doped Bi2212 compound. 49 The reduction in FWHM 
on going away from the node is observed along both 
the bonding and antibonding quasi-particle bands. This 
experimental momentum dependence is consistent with 
the off-diagonal results shown in Fig. 7 (a), but is diffi- 
cult to be explained by the magnetic excitation around 
q = (ir, 7r) and the diagonal electron-phonon coupling. 
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FIG. 6. Temperature dependence of the imaginary part 
of the memory function ImM(u>) at u> = (bold line). 
Broken line indicates a curve of ImM(0) = 2T. The di- 
agonal and off-diagonal electron-phonon coupling constants, 
(doff, Hdia), are chosen to be (0.03eV, 0). Other parame- 
ters are lu p h = 0.07eV, Ao=0 and t = 0.3eV. We use that 
Xo = 0.81ta 2 calculated for the half-filled tight-binding model 
in the two-dimensional square lattice with NN electron hop- 
ping. 
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FIG. 7. Contour map of the real part of the fermion renor- 
malization factor ReZ(k,Lj — Ao) in the off-diagonal coupling 
case (a), and in the diagonal coupling case (b). Broken lines 
indicate the fermi surface in £r. The diagonal and off-diagonal 
electron-phonon coupling constants, (g ff, Hdia), are chosen 
to be (0.03eV, 0) in (a), and (0, 0.02eV) in (b). Other pa- 
rameters are ui p h = 0.07eV and A = 0.035eV. 
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III. VERTEX CORRECTION 

A. Formulation of the Vertex Correction 

We now go beyond the slave-boson mean field treat- 
ment in the electron-phonon interaction in t — J model. 
This is carried out by taking into account the fluctuations 
around the mean fields, leading to the vertex correction 
in the fermion-phonon interaction. We formulate this ef- 
fective electron-phonon coupling in the SU(2) slave boson 
theory 35-37 which imposes that the fermion operators 



fil 

-A 



(28) 



are taken to be the local SU(2) doublet. This is an ex- 
act symmetry in the t — J model at half filling 35 where 
the two order parameters, i.e., the fermion hopping Xij( = 

E CT fi*fia) and thc fermion paring Ay(= fyfji-fufji), 
become equivalent. Even at the finite but small dop- 
ing, it is expected that the low energy physics is cap- 
tured by this picture. The generalized SU(2) theory 
away from the half filling is provided by introducing 
thc SU(2) doublet for bosons hi = {h li ,h 2i ) t . 36 Then, 
the physical electron operator is Ci G — h^i^ with con- 
straints \ ^la^ia +h\fhi = 0. Three bosonic fields 
a\(l = 1 ~ 3) are required as Lagrange multipliers to 
keep these constraints. It is known that the Lagrangian 
in the SU(2) theory is rewritten as being similar to that 
in the U(l) theory. This is carried out by representing 
the two-component bosons as a SU(2) rotation, denoted 
by the rotation matrix gi, of the standard slave-boson 
hi. 38. By integrating out the variable gi, the Lagrangian 
for the t — J model in the SU(2) theory associated with 
the electron-phonon interaction is obtained as 



L — L t j + L e i- p h, 
where L t j is the electronic part: 

to'} (*j),<r 



(29) 



LtJ = ^E Tr[C/* Uij] + J - E (^lUij^ja + H.c) 



E h\ (d T - ia 3 + n) hi 

i 

t E (flfi'W) + H - c ) > 



to),. 

with a 2x2 matrix 



Uij = 



Xiji 
A*-, 



Xij 



(30) 



(31) 



The second term in Eq. (29) is the the fermion-phonon 
interaction: 



k,q,f-i,cr 

xl(b%+bf), (32) 



which is nothing but thc Hamiltonian Eq. (15). From 
now on, we focus on the off-diagonal electron-phonon 
coupling in g^{k,q) of the present interest. 

The fluctuations around the mean-field saddle-points 
in the SU(2) theory are examined, in detail in Ref. 37, as 
the low-laying collective modes in the superconducting 
ground state. For the fermion hopping and the fermion 
paring, as well as the slave boson and the Lagrange mul- 
tipliers, we introduce small deviations from the saddle- 
points; 



Xij = Xo + Sxtj, 
Ay = A (-l)r» + <SA y) 
hi = r (l + Sri), 
ia\ = 6 l3 \ + i8a\. 



(33) 



rji = (1,-1) for a direction of the bond I = (x,y). We 
adopt the radial gauge where hi is chosen to be real and 
its phase degree of freedom is absorbed in the Lagrange 
multipliers. Up to the quadratic fluctuations, we have a 
fermionic part of the Lagrangian as 



J MF 



Li 



(34) 



with the constant term L . Lmf is the mean- field part 

Lmf = ~27v^I *L (^ To - ^ T3 - A £ T1 ) ( 35 ) 



where we define the mean-field quasiparticle energy 



(36) 



with X£ = 2Jxo(cosfc a; +cosfc2 / ) and A^ = 2 JA (cos fen- 
ces fcy). The conventional definition of the supercon- 
ducting gap A corresponds to 2JA . The bare fermion 
energy is considered up to the next NN hopping as 
tr = 2rQ{t(cosak x + cos ak y ) + t\ cos ak x cos ak y } . The 
mean- field solutions, xo, Ao and Ao, are determined by 
the saddle-point equations: 



k k 

k k 



(37) 
(38) 
(39) 



with 



/?£ = cosk x — cos ky and 7^ 



cosk x + cosk v . The calculations in Sec. Ill are based 
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FIG. 8. Feynman diagram for the effective electron-phonon 
interaction vertex. Bold lines, double lines and broken lines 
indicate the fermioin propagators Go(k, iu) n ), the bare collec- 
tive-mode propagators Gjo(q) = Mj(<j) -1 , and the phonon 
propagators D(q,iu> m ), respectively. 



on Lmf + L e i- p h where £g is obtained by the tight-bind 
fitting in spite of Eq. (36). The linear coupling of the fluc- 
tuation with the fermionic degree of freedom is described 
by the second term in Eq. (34): 



Li = 



'2N i _ 

k,q,a 



^l+qa fk,k+<r) ^ ke- 



rn 



(Xq) v and (f%% + ~)v (f = 1 ~ 12) arc the fluctuations 
with 12 components and their structure factors, respec- 
tively, defined as 



(*4x> <v SA '^ SA %' 5A '^ 5A '?v 
S x'g X , Sx'qy, 5a\, Sal, Sal, Sr^j, 



(41) 



and 
Trr - 

J k,k+q 



^ — 2JS X T a , —2JSyT , 2JC X T\, 2JCyTl, 

—2JC x t 2 , —2JC y r 2 , —2JC x t 3 , —2JC y T3, 



-ITi, -IT 2 , ~IT 3 . 



-(*s + W T:i 



(42) 



with abbreviations Si = sin(fc; + qi/2) and Ci — cos(fc; + 
qi/2) for I = (x,y). S X qi(= S X ' ?l + iSx'^) and SA^(= 
SA'^+iSA'^) are the Fourier transforms of S\ij and SAij, 
respectively, where I indicates a direction of the bond 
connecting site i and site j. The third term L 2 in Eq. (34) 
is the quadratic term of the fluctuations. Explicitly, 

Li = -2irtJ2 Sa l Sr -k 

k 

- r oX! (^fe^ 1 +^iCosa/c :E cosafc y A 2 - A ) \5r^\ 2 



+ ^_E (l*xal 9 + l«A a | a ), 

k,l=(x,y) 



(43) 



with Ai = - Eg and A 2 = -2 £r 

The effective fermion-phonon coupling is obtained in 
this scheme. The dynamics of the fluctuations x$ is gov- 
erned through the interaction between the fluctuations 
and fermion. By integrating out the fermionic degree 
of freedom in the action, we obtain the RPA-type sum- 
mation of the bare Green's function for the fluctuations 



£r: cos ak x cos ak y 
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FIG. 9. Contour map of (a) the effective paring interaction 
F(q) in the off-diagonal coupling case with the vertex correc- 
tion, and (b) the vertex function F(q, 0). Parameters are cho- 
sen to be t/J = 3, ti/J — —1.5, Hoff/J = 0.1, u>ph/J = 0.5 
and x — 0.1. Numerical data are plotted in a unit of J. 

Gja(q) = F.T.(xg-(T)x_JO))o, where (• • -)o indicates an 
expectation in L 2 , and the fermion polarization function 
Mp(q,iuJn) defined by 



M F {q,iuj n ) = .4+ g -,fe G o(Mw„ 



X G (fc + q, iuJ m + iVrOfk+frk- 



(44) 



Go(k,iuj m ) is the bare fermion Green's function calcu- 
lated by Lmf (Eq. (35)). The effective fermion-phonon 
vertex is now given by the bare interaction plus the cor- 
rection due to the collective modes of the fluctuations 
(Fig. 8). Explicitly, the coupling constant is 



g^(k, q; iw n ) = rft ■ A(q, iuj n ) ■ /g £+f . 
ffi is a vector with 12 components defined as 



(45) 



2NMuj i i 



: sin 



(I) 



{^ 7 cos(|)+cV 8 cos(|)}. (46) 
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FIG. 10. Superconducting condensation energy AE. Solid 
and broken lines indicate AEui p h/ J 2 with and without the 
vertex correction, respectively. The inset shows the real part 
of the fermion renormalization factor KoZ(k, 0). Parameters 
are chosen to be t/J = 3, ti/J = —1.5, u) ph /J — 0.5 and 
x = 0.1. 



The finite components at I = 7 and 8 in rft are attributed 
to the fact that the fermion-phonon vertex is proportional 
to T3. A(q 7 ioj n ) is the energy-dependent vertex function 
of a 12 x 12 matrix 



A(q,iu> n ) = 



Mj{q) 



Mj(q) + M F (q, iui n ) ' 



(47) 



represented by the inverse of the bare fluctuation Green's 
function Mj(q) = Gjo(q)^ 1 and the rcnormalized fluc- 
tuation Green's function 



Gj(q,iuj n ) 



Mj{q)+M F (q,iu n y 



(48) 



which describes dynamics of the collective modes in the 
superconducting ground state. 



B. numerical results with vertex correction 

We show, in this section, the numerical results ob- 
tained by taking into account the vertex correction for the 
electron-phonon interaction in Eq. (45). As shown below, 
the effects of the vertex correction depends strongly on 
the physical quantities. Some of them are seriously mod- 
ified, while others are insensitive to the vertex correction. 
The parameter values are given in a unit of J, and we 
use t/J = 3, ti/J = —1.5, g ff/J = 0.1, x = 0.1 and 
^ph/J = 0.5. The holon is assumed to be condensed 
and its amplitude ro is taken to be y/x. The obtained 
mean-filed solutions are A o =0.278 and \ = 0.370. 

The vertex correction is a crucial for the calculation 
of the paring function and the condensation energy. The 
momentum dependence of the effective paring interaction 
F(q) is obtained by Eq. (20) where g^(k,q) is replaced 



by the effective coupling constant 7j^(k,q;iui n — 0) de- 
fined in Eq. (45). The numerical results are presented in 
Fig. 9 (a). The vertex correction does not only enhance 
magnitude of the attractive interaction, as will be dis- 
cussed in more detail below, but also shifts the minimum 
region from (tt, 0) to (tt(1 - <5),0) with 6 = 0.3. This is 
exactly the momentum region where the minimum of the 
anomalous dispersion is observed in the inelastic neutron 
scattering. 14 ' 15 We plot, in Fig. 9 (b) the contour map 
of the absolute value of the vertex function T(q,u) = 0) 
defined by 

T(q, iuj n ) = A(q, iw n ) 7 „ cos ) 

+A(ftiw n ) 8 „cos(|^}. (49) 

We resolve this vertex function into the contribution from 
the collective modes by following the classification in 
Ref. 39. The strong intensity around q = (0, 0) in T(q, 0) 
is attributed to the Goldstone mode, i.e. the uniform 
collective phase mode of the superconducting order pa- 
rameter. This eigen vector is denoted as 



Sauq 



6A gy> 



(50) 



where a indicates the overall phase of \ij and A^- around 
the staggered- flux order parameters. This contribution 
in the effective fermion-phonon interaction is, however, 
irrelevant, because of the factor sin(g M /2) in Eq. (46) and 
also because the long range Coulomb interaction lifts up 
the frequency to plasma energy in real system. A large 
contribution around (it, 0) is mainly attributed to the two 
collective modes, i.e. the amplitude mode 



6Am = AoSA'^ + xoSx'qi, 



(51) 



for I = (x, y) representing the amplitude fluctuation for 
the order parameter \xqi\ 2 + |A«j-;| 2 , and the transverse 
and longitudinal <f) gauge modes 



hq = XoSA' ?y - IA 5x'qy, 



(52) 
(53) 



The amplitude fluctuation SAqi has a large intensity 
around (tt,0) and (0,7r), and increases the effective par- 
ing interaction F(q) around these momenta. This may be 
a remnant of the instability toward the eolumnar-dimer 
state which is well known in the insulating case. 50 ' 51 The 
4> modes are the amplitude modes with the out-of-phase 
fluctuation for the fermion hopping Xij and paring Ay. 
The longitudinal <\> gauge mode tpLq shows a large in- 
tensity in the low energy region around q = (tt/2,0), in 
addition to the peak structure around q = (ir/ '2, n) sug- 
gested in Ref. 39. This mode shifts the minimum point of 
F(q) from (tt, 0) without vertex correction to (tt(1 — <5), 0) 
with 5 = 0.3 in the low energy region. 
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FIG. 11. Real part of the phonon self-energy without 
the vertex correction (a), and with the vertex correction 
(b). Parameters are chosen to be t/J = 3, ti/J = —1.5, 
Hoff/J = 0.1, LOph/J = 0.5 and x = 0.1. 



The superconducting condensation energy AE is given 
by the expectation value of the effective BCS Hamilto- 
nian, i.e. Eq. (20) where the vertex correction is included 
in the coupling constant. The numerical results of AE as 
function of the bare off-diagonal coupling constant g a ff 
is presented in Fig. 10 associated with the renormaliza- 
tion factor ReZ(k, u> — 0). As shown in the broken line 
in Fig. 10, the condensation energy is positive in the case 
without vertex correction. Thus, the ci-wave supercon- 
ductivity due to the present LO phonon is stabilized by 
introducing the vertex correction. By considering the ve- 
locity ratio R v across the kink energy, a value of g f / / J 
is expected to be 0.2 where KcZ(q, 0) is about 0.6 <~ 0.8 
at q = (0,0) and (7r,7r), and AEuj ph /J 2 = -0.02 is con- 
verted into about 40K, taking J = lOOmeV. The actual 
condensation energy is of the order of 7K for optimally 
doped YBCO 52 , which is smaller than this value. One 
can consider various reasons for this discrepancy, (i) the 
estimation of the condensation energy has an uncertainty 
due to the subtle subtraction of the phonon contribution, 
(ii) there are other mechanisms to reduce the pairing and 
hence the condensation energy which are not considered 
here, and (hi) the determination of KeZ(k,0), i.e. the 
electron-phonon coupling constant A from the experimen- 
tal velocity ratio may be overestimated. 

The momentum dependence of the phonon self-energy 
is presented in Fig. 11. The phonon self-energy 
Y, p h(q,u) = uiph) is calculated from the fermion bubble 
of the bare fermion Green's function Go(k, u>) associated 
with the vertex correction. We consider the phonon mode 
where the oxygen ions vibrate along the x direction. As 
shown in the Fig. 11(a), the softening of the phonon dis- 
persion is remarkably seen along the (q x , 0) direction with 
1y ^5 7r/2. This is the region where the softening of the 
oxygen vibrational LO phonon mode is reported in the 
inelastic neutron scattering experiments. The phonon 
softening is gradually reduced with increasing q y , and 
along the {q x ,^) direction, the phonon frequency is al- 
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FIG. 12. Contour map of the fermion spectral density 
A(k, u>) along the nodal direction without the vertex correc- 
tion (a) and (b), and with the correction (c) and (d). (a) and 
(c) are for the spectra with Ao = 0, and (b) and (d) are for 
the spectra with A/ J = 2A /J = 0.57. Parameters are cho- 
sen to be t/J = 3, ti/J = —1.5, Hoff/J = 0.1, u>ph/J = 0.5 
and x — 0.1. The dispersion relation of the quasi-particle 
energy £g is chosen to be the same form with that in Fig. 3. 
Numerical data are plotted in a unit of J -1 . 



most unrenormalizcd. By introducing the vertex correc- 
tion (Fig. 11(b)), these tendency become significant. Al- 
though the qualitative feature of the calculated phonon 
softening is consistent with the experiments, absolute 
value of the frequency shift is much larger than the ob- 
servation, in particular with the vertex correction. This 
may indicate that the adopted electron-phonon coupling 
constant is overestimated, as mentioned previously, and 
a further refinement for the vertex correction is required 
for more quantitative discussion. 

On the other hand, the electronic spectra are not sen- 
sitive to the vertex correction as shown below. Figures 12 
show the ARPES spectra along the nodal direction with 
and without the vertex correction. The energy depen- 
dence of the vertex function is taken into account in the 
calculation of the spectra. Comparing Figs. 12 (a), (b) 
and (c), (d), it is concluded that the vertex correction 
does not play significant roles in the ARPES spectra. 
This is because the electronic self-energy is given by the 
contributions from the wide range of the phonon momen- 
tum q, and the vertex correction is large only in a rather 
limited q region as shown in Fig. 9(b). This fact is also 
consistent with the shift of the kink position in the su- 
perconducting state; the kink energy shifts from ui p h to 
uiph + 2aA with a = 0.7. However, this does not repro- 
duce the experimental observation that the kink position 
does not change across the superconducting transition. 

We present, in Fig. 13, a contour map for the renormal- 
ization factor for fermion KcZ(k,oj = 0) with the vertex 
correction. In the case where the vertex correction is 
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FIG. 13. Contour map of the real part of the fermion renor- 
malization factor ReZ(k, 0) in the off-diagonal coupling with 
the vertex correction. Parameters are chosen to be tj J = 3, 
ti/J — —1.5, ]}off/J = 0.1, ujph/J = 0.5 and x — 0.1. 

included, the essence of the momentum dependence of 
RcZ(fc, 0) pointed out in the previous section, i.e. the 
large renormalization around k — (0,0) and (7r,7r), re- 
mains, although ReZ(/c,0) around k — (7r,7r) is slightly 
reduced. 

IV. DISCUSSION AND CONCLUSIONS 

The interplay between the electron-phonon interaction 
and strong electronic correlation is one of the most impor- 
tant issues in many systems such as cuprates, mangan- 
ites, and Ceo- Especially in cuprates, it has been often 
claimed that the strong correlation dominates to hide 
or reduce the electron-phonon interaction. This seems 
to be true for the diagonal coupling of phonon because 
the charge density fluctuation is suppressed by the Mott 
physics. This is assumed also in the present paper ne- 
glecting the diagonal electron-phonon coupling. However 
it also plays an important role in the dynamics of the 
hole observed in the ARPES of the undoped cuprates 53 
because the photo-doped hole is subject to the rather 
strong coupling to the breathing phonon mode. We be- 
lieve that this electron-phonon interaction is the origin 
of the anomalously large width of the ARPES spectra. 54 

Now some remarks are in order on other works on 
the electron-phonon interaction in cuprates. Rosch and 
Gunnarsson 25 studied the electron-phonon interaction in 
the t — J model in a similar way estimating the value 
of the coupling constants from the r-dependence of the 
transfer integral and the inter-site Coulomb interaction 
between the O and Cu ions (r is the distance between 
these two ions). The off-diagonal coupling constant is 
an order of magnitude smaller than what we estimated, 
while the diagonal-coupling constant is huge of the or- 
der of 0.25eV. Therefore they concluded the diagonal 



coupling is the dominant mechanism of electron-phonon 
interaction in cuprates, even though it is reduced con- 
siderably by the same mechanism we discussed in this 
paper. However the small value of off-diagonal coupling 
constant is due to the rather subtle cancellation between 
two terms, and might change in a more realistic calcu- 
lation. Also the vertex correction strongly enhances the 
off-diagonal coupling. 

Huang et al. 23, recently studied the electron-phonon 
vertex correction for the diagonal coupling due to the 
Coulomb correlation, and found that the large momen- 
tum transfer scattering is reduced considerably while that 
of small momentum transfer is enhanced. This feature is 
consistent with the earlier analytic studies on the vertex 
correction. 40 The important point made in that paper 23 
is that the interaction between the quasi-particle and 
phonon is enhanced by the electron correlation because 
of the enhanced density of states. This might have an 
important implication for the kink structure observed in 
the (7r,0) momentum region. 34 

In summary, we have studied the electron-phonon in- 
teraction in HTSC cuprates motivated by the recent 
ARPES experiments. The off-diagonal electron-phonon 
coupling, which modulates the hopping integral of the 
Zhang-Rice singlet and the superexchange interaction, 
is relevant in the small doping region. We formulate 
the vertex correction beyond the mean-filed slave boson 
theory. As well as the slave boson and the Langrange 
multiplier, the fluctuations for the fermion hopping and 
paring lead to the collective modes in the superconduct- 
ing state. 39 Thus, the dynamics of the collective modes 
dominates the vertex function for the electron-phonon in- 
teraction. A characteristics in the off-diagonal coupling 
constant g^j(k,q) are consistent with several k (elec- 
tron momentum) and q (phonon momentum) dependent 
experiments; the anisotropy in the quasiparticle renor- 
malization reflecting the momentum-dependent velocity 
ratio across the kink and its peak-width, and the phonon 
momentum dominating the strong electron-phonon cou- 
pling. With the vertex correction, the attractive inter- 
action increases and the phonon momentum induces the 
interaction shifts to (tt(1 — <5),0) with S = 0.3. These 
changes are caused by the amplitude fluctuation mode 
and the (j> gauge mode, that is, the collective fluctua- 
tions of the order parameter amplitude \xij\ 2 + |Ay| 2 , 
and its out-of-phase fluctuation mode, respectively. How- 
ever, the kink position in the ARPES spectra shifts from 
ujph to ujph + 2aA across T c in our calculation, while it 
does not in the experiments. The resolution of this dis- 
crepancy is left for future studies. 

It is evident that the strong electron correlation plays 
the major role in the physics of HTSC cuprates, and the 
role of the phonon is also modified by it. Some specific 
channel of the electron-phonon interaction is selected to 
be enhanced by the strong correlation, and this will avoid 
the competition with the magnetic mechanism working in 
the different momentum regions. Therefore it seems that 
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the electron correlation and the electron-phonon interac- 
tion seems to collaborate to enhance the d-wave pairing. 
The role of the other phonons besides the half-breathing 
mode studied here, is left for future investigations. 
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